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PHOTOIONIZATION FEEDBACK IN LOW-MASS GALAXIES AT HIGH REDSHIFT 
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ABSTRACT 

The cosmic ultraviolet (UV) ionizing background impacts the formation of dwarf galaxies in the low-redshift 
universe (z 5a 3) by suppressing gas infall into galactic halos with circular velocities up to v c ; lc ~ 75kms _1 . Using 
a one-dimensional, spherically symmetric hydrodynamics code (Thoul & Weinberg 1995), we examine the effect 
of an ionizing background on low-mass galaxies forming at high redshifts (z J> 10). We find that the importance 
of photoionization feedback is greatly reduced, because (1) at high redshift, dwarf-galaxy sized objects can self- 
shield against the ionizing background, (2) collisional cooling processes at high redshift are more efficient, (3) 
the amplitude of the ionizing background at high redshift is lower, and (4) the ionizing radiation turns on when 
the perturbation that will become the dwarf galaxy has already grown to a substantial overdensity. We find that 



because of these reasons, gas can collect inside halos with circular velocities as low as v c 
This result has important implications for the reionization history of the universe. 



lOkms" 1 atz> 10. 



1. INTRODUCTION 

The recent detection of a large optical depth to electron scat- 
tering (j e = 0.17 ±0.04) by the Wilkinson Microwave Anisotropy 
Probe (WMAP, Bennett et al. 2003) suggests that the first sources 
of light significantly ionized the intergalactic medium (IGM) as 
early as redshift z ~ 17 ±4 (Kogut et al. 2003; Spergel et al. 
2003). A number of papers (see a recent review and summary 
in Haiman 2003) have tried to reconcile this observation with 
the preceding discovery of a significant global neutral hydro- 
gen fraction inferred from Gunn-Peterson troughs in quasars 
at redshift - 6.3-6.4 (Becker et al. 2001; Fan et al. 2003; 
White, Becker, Fan & Strauss. 2003). Haiman & Holder (2003) 
and Onken & Miralda-Escude (2003), in particular, have em- 
phasized the role that radiative feedback mechanisms can play 
in regulating the evolution of the global emissivity of ionizing 
radiation, and simultaneously accounting for these two appar- 
ently conflicting observational results. What is suggested by 
these observations is that reionization is completed at redshifts 
of z ~ 6-7, but that in addition, there is a significant "tail" 
of (perhaps only partial) ionization extending to much higher 
redshift (z ~ 17). 

If the reionization history is regulated by radiative feedback 
effects, such a tail of ionization can naturally be produced. The 
first ionizing sources form in low-mass dark matter (DM) halos 
with virial temperatures between ~ 10 2 and 10 4 K around z = 20, 
in which cooling by H2 molecules can bring the Jeans mass 
down to several hundred M Q . The UV photons produced by 
these early population of stars or miniquasars photo-ionize and 
photo-heat the gas surrounding them (Shapiro et al. 2003; Oh 
& Haiman 2003) and dissociate H2, the main coolant (Haiman, 
Rees & Loeb 1997). This limits the star formation rate and 
the contribution of these halos to reionization (Haiman, Abel & 
Rees 2000). 

The global star-formation can increase again once more mas- 
sive halos form, with virial temperatures of 7^ J> 10 4 K (or cir- 
cular velocities of v C ac <; lOkms" 1 ). In these halos, gas can 
condense and cool via an initial phase of atomic cooling (Oh 
& Haiman 2002). In practice, however, previous models of 
the reionization history have excluded star formation in ha- 



los up to much higher virial temperatures, with T vlr <; 10 5 K 
(vdrc Ji 50kms _1 ) in regions of the IGM that have already been 
ionized (e.g. Haiman & Loeb 1997, 1998; Wyithe & Loeb 
2003; Haiman & Holder 2003). The motivation for this ex- 
clusion is that the collapse and cooling of gas is thought to 
be significantly suppressed in these photo-heated regions. The 
critical halo size up to which such a suppression is important 
has been addressed by several authors (Efstathiou 1992; Thoul 
& Weinberg 1996, hereafter TW96; Navarro & Steinmetz 1997; 
Kitayama & Ikeuchi, hereafter KI00), who found a critical value 
of v c i rc ~ 50kms _1 . The exact value of this threshold plays a 
crucial role in the reionization history. This motivates us in the 
present paper to revisit the effects of the UV feedback at high 
redshift. 

It is important to realize that existing studies of the UV feed- 
back mainly focus on low-redshifts (z 3), addressing the sup- 
pression of the formation of dwarf galaxies in the relatively 
nearby universe. Whether or not the same feedback occurs at 
high redshift is especially important to clarify in light of the 
WMAP results. Specifically, a strong photoionization feedback 
can significantly delay the completion of reionization, provided 
the ionizing photon production efficiency is high in these halos 
(see Figure 4 in Haiman & Holder 2003). In general, one would 
expect photoionization feedback to have less of an effect at high 
redshift, because cooling times are shorter, the intensity of the 
ionizing background is lower, and halos are more compact and 
can be self-shielding. For example, KI00 included radiative 
transfer effects in one-dimensional simulations similar to those 
of TW96. They found that this lowered the critical velocities 
by ~ 5kms _1 , at low redshift, and that this decrease was more 
significant (~ 5kms _1 ) at high redshifts. In three-dimensional 
simulations, Gnedin (2000) found that the characteristic mass 
scale for suppressing the gas fraction in non-linear halos de- 
creases towards high redshifts, with a dependence that roughly 
corresponds to a fixed circular velocity of 20 km/s at all red- 
shifts (and provided a fitting formula that has been adopted in 
the reionization models of Cen 2003). 

In the present paper, we extend previous calculations of the 
photoionization feedback to high redshifts, using spherically 
symmetric simulations adapted from Thoul & Weinberg (1995). 
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Our goal is to quantify the suppression of gas infall into halos 
that form at z ~ 10, and to elucidate the physics that distin- 
guishes the effects of the UV background at z ~ 10 from that at 
z ~ 2. The recent WMAP results suggest that an ionizing back- 
ground is already in place at high redshift, and much of early 
structure formation takes place under the influence of an ion- 
izing background. Previous works nevertheless suggests a di- 
minishing importance of the background towards high redshifts 
(KI00, Gnedin 2000). In this paper, we confirm these results, 
and explicitly compute the change in the characteristic circu- 
lar velocity with redshift. Using a range of simulation runs, 
we also quantify the relative importance of different physical 
effects that determine the condensed gas fractions we obtain. 
Finally, we discuss the significance of (the lack of) feedback 
for the high redshift reionization history in light of the WMAP 
results. 

The rest of this paper is organized as follows. In § 2, we 
give a brief description of the code we adapted, together with a 
few modifications we made for the present application. In § 3, 
we describe the results of our runs, which are then discussed in 
more detail in § 4. In § 5, we summarize our conclusions and 
the implications of this work. Throughout this paper, we adopt 
the background cosmological parameters as measured by the 
WMAP experiment, O m = 0.27, fi A = 0.73, Sl b = 0.044, ^ = 0.71 
(Spergel et al. 2003). 

2. DESCRIPTION OF THE CODE 

The code we use is a modified version of a one-dimensional, 
spherically symmetric code (Thoul & Weinberg 1995, hereafter 
TW95). It evolves a mixture of dark matter and baryon fluids 
by moving concentric spherical shells of fixed mass in the radial 
direction. A description of the numerical method and details 
of the code can be found in the above reference. We start all 
runs with an initial density profile following the average shape 
of a 2cr peak in a Gaussian random field (Bardeen et al. 1986; 
BBKS). The initial density profile is sampled by 6000 dark mat- 
ter and 1000 gas shells out to the outermost shell at r « 2.3r/, 
where jy is the filter radius for the BBKS profile. We have ex- 
plicitly verified that increasing and decreasing the number of 
shells by a factor of three changes our answers for the cold gas 
fraction below by at most 5 percent. This will introduce an 
error in the circular velocities of at most ~ 3 km/s. The total 
mass (baryons + dark matter) of the object collapsing at red- 
shift z = 1 1 is denoted by Mf. This mass is initially contained 
in a sphere of radius Iff, in which we fix the mean overdensity 
to be 5 = (p-p)/p = 0.136. Note that this places the redshift of 
collapse at z c ~ (0.136/1.69) x (1+z,) — 1, where zi is the ini- 
tial redshift. We chose z, based on the desired collapse redshift 
(e.g., Zi = 160 for z c = 1 1). The total mass of the collapsing ob- 
ject in our simulation is 1 .5Mf (see TW95 and TW96 for more 
details). 

The code was originally written for a standard cold dark mat- 
ter (SCDM) cosmology (fl m = 1 , A = 0, H = 100 km/s/Mpc). 
In this paper, we adopted the ACDM cosmology favored by 
the WMAP results. However, we are concerned with the col- 
lapse of objects at high redshift, where the dynamical effect of 
the dark energy is negligible. In practice, we therefore needed 
to change only the ratio n m /fij (which sets the number ratio 
of dark matter and gas shells), and the time-redshift relation 
(f = 425[(1 +z)/12]" 3 / 2 Myr in the redshift range 2 < z < 160). 
Because the time-scale for Compton cooling off the cosmic mi- 
crowave background (CMB) is shorter than the Hubble time for 



z ^ 6.5, we added Compton cooling to the other cooling terms 
in the code (e.g., Katz, Weinberg, & Hernquist, 1996), 

A C ompt= 1-017 x 10- 37 (r c -r cmb )r c 4 mb « c erg sec" 1 cm" 3 , (1) 

where n e is the number density of electrons, T e the tempera- 
ture of the electrons (equal to the gas temperature), and T cm b = 
2.725(1 +z) is the temperature of the CMB (Fixsen & Mather 
2002). Finally, it will be useful to recall the relation between 
the circular velocity v c ; rc of a dark matter halo and its virial 
temperature: 

iw - 3a xio-(A) (2) 

Below, we study the collapse of galaxies with circular veloci- 
ties in the range 10-100 km/s, for which the mean molecular 
weight is fi = 0.59. Note that below ~ 10 km/s, H 2 molecule 
formation and cooling has to be considered. The effect of the 
UV background on such minihalos has recently been studied 
elsewhere (Oh & Haiman 2003). The circular velocity of halos 
is related to their mass and collapse redshift (in the absence of 
any pressure or radiation field) as 

v - =18A {wm-) nrj (aos ) km/s ' (3) 

where Mf is the mass of the halo that collapses at redshift z c . 

3. RUNS AND RESULTS 

Although the WMAP data suggest that there is an ionizing 
UV background at z ~ 17, its evolution is poorly constrained. 
For simplicity, in our simulations we adopt a power-law ion- 
izing background that is described by two parameters, /„ = 
h\(ylvi)~ a x 10" 21 erg/s/Hz/cm 2 /sr, where v L is the Lyman 
limit frequency. In addition, the postulated ionizing background 
is assumed to turn on suddenly at the location of the collapsing 
object at a redshift z on . This sudden turn-on is appropriate be- 
fore reionization is complete, i.e. when the Stromgren surface 
from a single ionizing source passes by a newly forming halo 
(see more discussion of the adopted background below). We 
consider two different amplitudes (/21 = 1 or/21 = 10" 2 , justified 
below) and the spectral slope a = 1, characteristic of quasars 
or massive, metal-free stars. The choice of a = 5 would be 
more appropriate for a normal stellar population (TW96), but 
since making the spectrum softer would only reduce the impor- 
tance of feedback further, we are conservative in our choice of 
a = 1. For the turn-on redshift, we adopt z on = 17 (motivated 
by WMAP ) or z on = 00 (to obtain a maximal effect). We run 
each simulation to redshift ~ 6. In each simulation run, follow- 
ing TW96, we define the fraction / coU = M c /M c (p = 0) of the 
baryonic mass that has collapsed (M c ) at redshift z = 11, rela- 
tive to the amount of baryonic mass that would have collapsed 
in absence of pressure (M c (p = 0) = (Clb/Cl m )Mf). Since we are 
considering halos well above the cosmological Jeans mass, in 
the absence of a radiation field we expect / co u = 1. A back- 
ground flux has an increasingly larger effect on smaller halos, 
so that in the presence of a radiation field, we expect / co n to 
increase monotonically with circular velocity. 

In Figure 1, we show the collapsed fraction as a function of 
circular velocity under various assumptions. To begin, we first 
reproduced the results of TW96, by adopting their redshift- 
evolution for the background flux, and considering halos col- 
lapsing at z c = 2 in their SCDM cosmology. Our results for this 
case are shown by the the circles in the upper left panel of Fig- 
ure 1 . The circles show that gas infall in systems with circular 
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velocities of ~ 35 km/s is suppressed almost completely, while 
the collapse of structures with circular velocities of 70 km/s is 
suppressed by 50% (Note that the often quoted value of 50 km/s 
is obtained by defining f co u at the time 2t c . Here we consider 
lf e instead, as explained in § 4.3). 

Hereafter, we will denote these two velocities by vo and V1/2, 
respectively. To facilitate comparison with previous studies, 
throughout this paper we focus on V1/2. which conveys the size 
of a halo that is strongly affected by the radiation. In the context 
of reionization, a question that may be more relevant is whether 
any gas may be available for star formation in the halo. The 
answer to the latter question is better conveyed by vo, for which 
we also quote values below. 

The rest of Figure 1 describes collapse calculations in the 
ACDM cosmology. First, the triangles in the upper left panel 
show the main result of this paper, and correspond to halos that 
collapse at z c = 11. We use the set of parameters J21 = 10" 2 , 
a = 1, and z on = 17. As the figure clearly shows, the values of vo 
and V1/2 change dramatically relative to the low-redshift case: 
we obtain vo « 15 km/s, and V1/2 ~ 20 km/s. In the remaining 
three panels of Figure 1, we investigate the physical reasons for 
the change in the critical circular velocities, by varying each 
parameter in succession. Our approach is to start from a fiducial 
set of parameters that describe the typical low-redshift case: 
z c = 2, J21 = 1, a = 1, and z on = 00 (in practice, z on = z ; ), and 
then change each parameter in succession, until we arrive at the 
high redshift model described by z c =11, ^21 = 10 2 , andz on =17. 
The order in which we change the parameters is arbitrary (and 
could somewhat modify the relative importance of each effect). 
Gas Density: Cooling through two-body processes, such as re- 
combination and collisional excitation cooling, is more efficient 
at high densities. The influence of increasing the density by a 
factor of (12/3) 3 = 64 on the collapsed fraction is shown in the 
upper right panel, by adopting z c = 1 1 (stars) vs z c = 2 (empty 
circles). The other parameters are held at their fiducial values. 
As the panel reveals, the density increase shifts the entire / co u 
curve by <~ 10 km/s to the left. In particular, v\n is decreased 
from ~ 80 to ~ 65 km/s (a factor of ~ 1 .2). 
Radiation Intensity: At redshift z J> 1 1, the UV background is 
expected to be significantly lower than at redshift z <~ 2. We 
therefore next decrease the amplitude of the flux by a factor of 
100 to /21 = 10~ 2 . This value is motivated below in § 4.1. The 
lower left panel shows the result of this decrease (asterisks) vs. 
721 = 1 (stars, reproduced from the upper right panel). Lowering 
the flux yields an additional factor of <~ 1.4 decrease in V\i 2 
from ~ 65 to <~ 45 km/s. 

Turn-on time of ionizing flux: A primordial perturbation that 
collapses at redshift z ~ 11 is already turning around at z ~ 18. 
As a result, one expects that such a perturbation evolves sig- 
nificantly without the presence of any UV flux before any flux 
turns-on (see discussion below). To include this effect in a sim- 
ple way, in the lower right panel we change z on = 00 (empty cir- 
cles, reproduced from the lower left panel) to z on = 17 (triangles; 
we have arrived back at the "original" curve shown in triangles 
in the upper left panel). The late turn-on results in an additional 
factor of <~ 2.2 decrease in V\n from ~ 45 to ~ 20 km/s. Note 
that photoheating of halos that collapse at redshift z~2 (such 
as the cases considered in TW96) starts well before they turn 
around, and the precise turn-on redshift is unimportant. The 
strong influence of the turn-on time in our high-redshift col- 
lapse case can also be seen directly in the evolution of shells in 
the upper pair of panels in Figure 2. 



In Figure 2, we investigate the various physical effects de- 
termining the strength of the UV feedback in more detail. The 
figure shows (upper to lower panels, respectively) the evolution 
of individual shells, together with their temperatures, ionization 
fractions, and cooling times, during the collapse of a halo with 
Vdrc = 32 km/s. The set of panels on the left assume that the UV 
background is switched on at z on =17, and the panels on the 
right describe runs in which the flux is on during the entire sim- 
ulation. All quantities are shown as a function of cosmic time 
normalized to the collapse time t c = 425Myr (the age of the uni- 
verse at z = 11; the upper horizontal axes show corresponding 
redshifts). We next describe the individual panels. 
Evolution of gas shells. The uppermost row of panels shows 
the evolution of two gas shells, initially enclosing masses of 
0.5M/ and Mf (we will hereafter refer to these shells as shell 
0.5 and shell 1). The radii are shown in units of the filter ra- 
dius Yf. The pair of solid curves shows the evolution of shells 
in the absence of any pressure in both panels. Shell 1 collapses 
approximately at t/t c = 1 as desired. 1 The dotted and dashed 
curves in the right panel correspond to /21 = 10~ 2 and /21 = 1. 
respectively, showing that the shells collapse later (or do not 
collapse at all) as the the intensity of the UV background is 
increased. In the left panel, the thick dotted line shows the evo- 
lution of the shells when the UV background (with /21 = 10" 2 ) 
is turned on at redshift z = 17. At this redshift, both shell 1 and 
0.5 have already turned around. The figure shows that shell 0.5 
was further along its evolution that shell 1, and therefore the 
background has a smaller influence in its subsequent evolution 
than in that of shell 1 . 

Evolution of gas temperature. The second row of panels from 
the top shows the evolution of the temperature of gas shell 1, 
with the dotted curve in the left panel representing the cases 
with /21 = 10~ 2 and the dashed curve on the right when J 2 \ = 1. 
For reference, we also show the CMB temperature (dot-dashed 
curves), and the adiabatic temperature evolution (oc (1+z) 2 , 
solid curves). Initially, heating takes place through adiabatic 
compression of the gas shells. Some heating will take place 
through inverse Compton scattering of the CMB photons by 
electrons, however this is very inefficient. We set the fraction 
of free electrons to n e /{nm + «hii) = 10~ 4 (the approximately 
constant value expected after decoupling, e.g., Peebles 1993). 
The left panel shows that the gas temperature is locked between 
the adiabatic the CMB temperature and an adiabatic evolution, 
until the background turns on. At z on , the shell is rapidly photo- 
heated to 10 4 K, irrespective of the amplitude of /21 (this is en- 
forced by the assumption of ionization equilibrium in our code; 
however, this assumption is justified since the photoionization 
time-scale is short, ~ 10 6 yr for 7 2 i = 10" 2 ). 
Evolution of ionization fractions. The third row of panels from 
the top shows the evolution of the neutral (solid curves) and ion- 
ized (dashed curves) fractions of hydrogen for shell 1 . The left 
panel shows the model with /21 = 10~ 2 and with the UV back- 
ground switched on at z = 17. As mentioned above, the free- 
electron fraction is set to 10~ 4 before the UV background turns 
on. Once the flux turns on, the gas is highly ionized. However, 
as shown by the solid curve in the left panel, significant neutral 
hydrogen remains in this case. This implies that collisional ex- 
citation cooling can dominate over other cooling processes, and 

'The exact collapse time is sensitive to the initial velocities assigned to the 
shells. We adopted velocities derived from linear theory for the growing modes: 

Vj = ffr,(l - y ), where (5, is the mean initial overdensity enclosed by shell i (see 
TW96 for details). 
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also that self-shielding effects can be important (see discussion 
below and in § 4.2). For clarity, the right panel only shows the 
case with /21 = 1, and demonstrates that the halo gas is quickly 
highly ionized, with the neutral fraction declining due to the 
dilution of the density in the expanding shell. 
Evolution of cooling timescales. The dominant cooling mecha- 
nisms are Compton, recombination and excitation cooling. The 
timescales for all processes are shown in the bottom row of 
panels for shell 1. The left panel shows that for shells that 
do collapse, collisional excitation cooling (mostly by Lyman a 
emission) dominates over recombination and Compton cooling. 
This is because the densities of the shells shown in the left panel 
are higher (the UV flux is lower), and the fraction of neutral hy- 
drogen is higher, than in the corresponding curves in the right 
panel. The right panel shows that for the /21 = 1 case, Compton 
cooling dominates at high redshift, followed by recombination. 
Collisional cooling of neutral hydrogen is so inefficient that it 
is off the plot. The fact that Compton cooling is dominating 
over the other cooling mechanisms is because this shell never 
collapses and remains at relatively low densities. In this case, 
Bremsstrahlung also contributes to the total cooling rate at the 
percent level, as shown by the dot-dashed curve in the lower 
right panel. 

4. DISCUSSION 

The main result of the previous section is that the suppression 
of gas infall by photoionization heating is much less severe in 
halos collapsing at z = 11 than at z = 2. At redshift z = 2 the 
presence of a UV background photo-heats the baryons to ~ 10 4 
K, which prevents a fraction of the baryons from condensing 
inside dark matter halos. This fraction is 50% for halos with a 
circular velocity of V\/ 2 ~ 75 km/s, in agreement with previous 
results by TW96. 

We found that for halos collapsing at redshift z = 1 1, the crit- 
ical velocity is reduced to V\/2 = 20 km/s. This is a change 
in V1/2 by a factor of 3.5, corresponding to a decrease in the 
halo mass by a factor of 3.5 3 « 42 (eq. 3). We found that sev- 
eral factors contribute to this reduction. The largest factor of 
« 2.2 came from our assumption that for halos collapsing at 
Z = 11, the UV background is not always "on" - we choose to 
turn on the flux when the perturbation has just turned around 
(see the lower right panel in Figure 1). This corresponds to 
Zon =17, the redshift where WMAP indicates that a significant 
fraction of the volume of the IGM is ionized. At high redshift, 
the densities are higher, and therefore cooling is more efficient. 
This was responsible for a factor of ~ 1.2 (see the upper right 
panel in Figure 1)). A comparable factor (« 1.4) came from 
lowering the flux /21 from 1 to 10" 2 (in the next section, we 
discuss our choice for the amplitude /21 = 10" 2 ). In compar- 
ison, TW96 found that the influence of photoionization came 
from the suppression of turnaround and collapse by heating, 
rather than from reduction in cooling rates, and that therefore 
the amplitude of the background was unimportant (though the 
slope was significant because of its influence on temperature). 
This agrees well with our conclusion that the late turn-on of the 
background flux (relative to f c ), accounts for most of the reduc- 
tion in V1/2. However, in our high redshift case, the effect of 
photoionization on cooling does have an additional impact, and 
the amplitude of the flux is therefore important, as well. 

Two of the simplifying assumptions in our treatment are that 
the gas is optically thin and in ionization equilibrium. In real- 
ity, when the collapsing gas is illuminated by a flux that turns 



on suddenly, the initial heating rate can significantly exceed the 
values we assume. This is because neither the extra heat input 
from the initial ionization of each atom (e.g., Hui & Gnedin 
1997), nor the temporary hardening of the spectrum due to ab- 
sorption by the large initial column of neutral atoms (Abel & 
Haehnelt 1999), are captured by our treatment. However, ef- 
ficient cooling at high redshift tends to mitigate these effects. 
In addition, we note that the gas in our high-redshift collapse 
runs, is never highly ionized, and so the equilibrium heating 
rates are not significant underestimates. In order to explicitly 
verify these assertions, we performed a run analogous to the 
fiducial high-redshift run shown in the left panels in Figure 2, 
except that ( 1 ) we artificially set the temperature at z on to 50,000 
K (the highest temperature to which the gas could get photo- 
heated due to self-shielding effects), and (2) then follow the 
non-equilibrium evolution of the ionization fractions and heat- 
ing rates. We find that these effects decrease / co n by less than 
1%. 

So far our main focus was to differentiate the value of V1/2 at 
z = 2 and z = 1 1 . In semi-analytic models, it is useful to know 
how V1/2 changes as a function of redshift. This is shown in 
Figure 3. The circles denote the case with /21 = 10~ 2 for all red- 
shifts (the error-bars show the range obtained using a different 
number of shells N DM = 2000,6000, 18000). Since this flux is 
not appropriate at lower redshifts, the triangles denote the case 
with the UV background evolving as a function of redshift The 
evolution for /21 is taken from KI00 (their eq. 4), and has /21 = 1 
for 3 < z < 6 and monotonically decreasing after this to ~ 10" 2 
at z <~ 14. This demonstrates that at lower redshifts z < 6 the 
value of v\ji moves up further, because of the higher value of 
/21 • Since in this evolution scenario, /21 is too high at redshifts 
5 < z < 6, the triangles may be viewed as upper limits. In § 4.2 
we will discuss that the the value of v 1/ / 2 will decrease, when 
self-shielding and radiative transfer are included. 

4. 1 . Ionizing Background 

As mentioned above, the evolution of the UV background 
is poorly constrained beyond z = 6. Under the assumption that 
the IGM is in photoionization equilibrium with a uniform back- 
ground flux, several authors have used the presence of the Gunn- 
Peterson trough in a redshift z = 6 .28 quasar (Becker et al. 200 1 ) 
to place the upper limit /21 0.03 on the background flux at 
z - 6 (McDonald & Miralda-Escude 2001; Cen & McDonald 
2002; Fan et al. 2002; Lidz et al. 2002). The background at 
higher redshifts would be expected to be lower; in general, un- 
der the assumption of photoionization equilibrium implies the 
relation 




between the neutral fraction (by number) of hydrogen atoms 
(1 -x), and the background flux J21. While 1 — x is believed to 
be ^ 10" 4 at z w 6, our assumption of /21 = 10~ 2 at z = 17 is 
equivalent to (1 -x) = 2.2 x 10" 3 at this redshift. 

A uniform background at z ~ 17, however, is likely to be an 
overly simplistic assumption. Indeed, our main motivation here 
for examining the photoionization feedback effect is to quan- 
tify the ability of halos with 10 kms" 1 <^ v c i rc <, 75 kms -1 to 
contribute to reionization - a question that would be less rele- 
vant once the discrete ionized regions percolate to establish a 
uniform background, implying that the IGM has already been 
fully ionized. As a result, a picture that is more relevant to 
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consider (and which is also more realistic), is that of a swiss- 
cheese topology for the IGM: discrete HII regions (or perhaps 
the overlapping cluster of a few HII regions) populate the other- 
wise still neutral IGM. In this case, we can estimate the UV flux 
seen by the gas inside an individual HII region. For example, 
let us consider a 10 8 M Q dark matter halo, in which a fraction 
/» = 0.1 of the baryons turns into stars, generating n p h t = 4000 
ionizing photons per stellar baryon at a steady rate over a pe- 
riod of t sys = 10 8 yrs, of which a fraction f esc = 0.1 escapes into 
the IGM (these fiducial numbers were adopted from Haiman & 
Holder 2003). At a distance r away from the central source, we 
then have the flux 

J (^)U)(tr) < 5 > 

f Q. b /Q. m \ / M \ / n phot \ / 10V\ 
V 0.16 ) V10 8 M Q / V4000/ V t sys )' 

where the fiducial radius r = 17.9 kpc is the maximum physical 
radius the HII region can obtain (in the absence of recombina- 
tion). More generally, the maximum physical radius the halo of 
mass M can ionize is 

0.022 n phot /. f esc M \ i/3 / 18 \ 

h h 2 40000.1 o.i io 8 m / \i+ z ) pc ' 

(6) 

It is most likely to locate another halo in the edge of an HII 
sphere, simply because the volume is largest there. We con- 
clude that the ionizing flux seen by a collapsing object at z ^ 10 
is small, J21 & 10" 2 , regardless of whether or not the individual 
HII regions have overlapped. 

It is also interesting to consider the effect of spatial correla- 
tions between ionizing sources and forming halos. The mean 
separation between halos of mass M=10 8 M atz= 17 is ~60 
physical kpc (Jenkins et al. 2001). Thus, the random chance 
of finding a second, similar-size halo within the HII sphere is 
~ 3%. However, the halos are highly correlated (with a bias pa- 
rameter b ~ 7 (Mo & White 1996), and a correlation function 
£(r = 20kpc) ~ 5), bringing the probability of finding a second 
halo within the HII region to ~ 10%. Furthermore, we are in- 
terested in would-be halos, in addition to fully formed ones. 
The mean separation between halos which are at their turn- 
around stage at z = 17 is ~ 3 times smaller than that of fully 
formed halos (~ 1. 5a peaks, rather than ~ 3cr peaks). One 
or two such halos in their formative stages should therefore be 
found in a typical HII region, experiencing photo-ionization 
feedback. This implies, in turn, that the fraction of halos po- 
tentially subject to the feedback is significantly larger than the 
simple "volume filling factor" of the HII regions. 

4.2. Optical Depth Effects 

Another important issue is that towards higher redshifts, ha- 
los with a fixed circular velocity can self-shield against the ion- 
izing radiation more effectively. If the gas is in hydrostatic equi- 
librium within an NFW halo (Makino, Sasaki, & Suto (1998), 
with concentration parameter c = 5), a line of sight across a halo 
with circular velocity v C i rc at redshift z has a total hydrogen col- 
umn density (from the center to the virial radius of the halo) 



^V H = 2.5xl0 21 (-^ 
V 15 kn 

fl b h 2 



circ 

kms -1 



0.135 



0.022 J V n m h 2 



1+z 
12 

1/2 



3/2 



cm 



(7) 



When such a halo resides in a photo-ionizing background, 
the gas is ionized and the column density of neutral atoms is 
reduced. In general, a simple expression analogous to equa- 
tion (7) does not exist in this case. Here we assume that the 
gas is optically thin (r = 0), and in ionization equilibrium with 
a background flux with amplitude J21. We then show the nu- 
merically computed optical depth (t = Nm^m, where (T H i is the 
photoionization cross section for hydrogen) for 15 km/s halos 
as a function of redshift in Figure 4. The solid/dotted curve 
corresponds to 72 1 = 1/1 — 2 (both at the temperature of 10 4 K). 
The curves show the optical depth from shell 1 inward to shell 
0.5, i.e. to the shell containing 50% of the gas mass, and are 
therefore relevant for computing the value of V\/2- The figure 
demonstrates that self-shielding can become important at high 
redshift: typical column densities at z > 12 are several order of 
magnitudes higher than at z ~ 2, both because of the increased 
densities and the lower amplitude of the ionizing background. 

Of course, NFW halos are collapsed structures, and therefore 
have densities > 178 times the mean density of the universe at 
any given redshift. Therefore they are more efficient self shield- 
ers than the objects in our simulations, which are exposed to a 
UV background at turnaround, when the mean density interior 
to shell 1 is only <~ 5 times the mean density. Furthermore, 
the curves in Figure 4 implicitly assume that within the virial 
radius of the dark matter halo, the mass fraction of baryons is 
equal to the universal value, i.e., that a fraction / co u = 1 of the 
baryons has collapsed in the dark matter halo. In our simula- 
tions runs, this fraction is reduced by the UV background to be 
less than unity, which further reduces the optical depth relative 
to the simple analytical cases shown in Figure 4. 

We therefore next measure optical depth directly in our sim- 
ulations runs. Since these runs do not include radiative transfer 
and assume the gas is optically thin, we are effectively check- 
ing whether these runs are self-consistent. We measure the 
optical depth in runs with z on =17, immediately after the UV 
background is turned on, and when the gas shells are just turn- 
ing around. We find that high-redshift halos can easily reach 
Nm 10 18 cm" 2 and can therefore self-shield even at the time 
of turn-around. More precisely, in our fiducial high-redshift 
runs with z c = H. J21 = 10" 2 , a = 1 and z on = 17, we find that 
at the redshift z on =17, the smallest halo we simulated (v c j rc = 
15 kms" 1 ), has a column density of 7V H i ~2x 10 18 cm" 2 between 
shell 1 (which roughly represents the edge of our object) and 
shell 0.5. The corresponding optical depth (t = 14) is shown 
as an open dashed circle in Figure 4. Likewise, the lower, 
solid circle in Figure 4 shows the optical depth at turnaround 
when /21 is changed to 1. As expected, the optical depth in 
both cases is lower than that for the gas in the NFW halo. 2 
As a result of the high optical depth (t = 14) in the 15 km/s 
fiducial high-redshift run, shell 0.5 will see only a fraction 
exp(-o-Hi/Vei) ~ 10" 6 of the UV background. This fraction is 
decreased by an additional factor of ~ 10 when the average 
transmission, < e~ T >= j- J dSle^ 8 '®, is considered. The 
above makes it clear that self-shielding is important at high red- 
shift, and the triangles shown in the upper left panel of Figure 1 
should therefore be regarded only as lower limits. 

2 When the optical depth to the center of the halos is calculated, a large con- 
tribution comes from the innermost shells, where the density, and therefore the 
neutral fraction, is high. However, the evolution of the innermost few shells is 
known to depend on the choice of the numerical boundary condition at r = 
(e.g., Haiman, Thoul & Loeb 1996; Forcada-Miro & White 1997), and we 
therefore chose to remain conservative, and to exclude these shells from our 
calculation of the optical depth. 
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To demonstrate the importance of shelf shielding more ac- 
curately, we performed one run which takes into account ra- 
diative transfer. We follow the approach of Kepner, Babul, & 
Spergel(1997). To keep the computing time relatively short, 
only the opacity due to hydrogen is included in a self consistent 
way (ignoring helium). For each shell, we calculate the trans- 
mission, averaged over solid angle, and obtain the processed 
spectrum and reduced photo-ionization and heating rates. Al- 
though this approach still underestimates the total opacity, it 
gives us much more stringent constraints on v x i 2 and vo- 

We show in Figure 1 how the point with v c ; rc =15 km/s shifts 
upwards from / co n = to 1, when the opacity due to hydrogen is 
taken into account in the simulations. Including radiative trans- 
fer further increases the optical depths of our halos, since the 
inner shells see (much) lower fluxes than we are assuming in 
our optically thin runs, and therefore collapse more rapidly and 
to higher densities. These conclusions are consistent with those 
of KI00, who apply solutions of the radiative transfer equation, 
obtained in a plane parallel geometry, to spherical 1-D simula- 
tions. This causes them to overestimate the true optical depth. 
They found that the value of V\/ 2 is decreased by 5 km/s at low 
redshift, and by as much as 10 km/s at higher redshift, as com- 
pared to the transparent case studied by us and by TW96. 

Finally, as argued above, the question of photoionization feed- 
back studied in this paper is most relevant to halos forming in- 
side existing HII regions of a single ionizing source. The UV 
radiation field is then very asymmetric. This is not, of course, 
captured in our spherically symmetric simulations. If, however, 
only one side of the collapsing halo is exposed to the UV radi- 
ation, this will make it easier for gas on the shadowed side to 
collapse. In summary, the issues of radiative transfer and three- 
dimensional topology are likely very important and should be 
addressed in future work. Here we simply note that our main 
conclusion - namely that feedback is less important towards 
higher redshifts - is likely to be significantly strengthened by 
both effects. 

Table 1 

Critical circular velocities (in km/s) for gas cooling 

DEFINED AT TWO EPOCHS: t c VS. 2t c AND IN DIFFERENT RUNS 
CORRESPONDING TO THE DIFFERENT PANELS IN FIGURE 1 (SEE 
§ 3 FOR DISCUSSION). 





\t c : 




2t c : 




change 


Vl/2 


vo 


Vl/2 


vo 


z = 2 data 


80 


45 


55 


30 


raising z (or p) 


65 


30 


37 


20 


lowering J2\ 


40 


20 


25 


15 


turn on at z = 17 


20 


15 


20 


15 


self-shielding 


10 


10 


10 


10 



4. 3 . Clustering and Merging 

Another issue is the epoch at which we define to look at how 
much gas has cooled and collapsed in a halo. An often quoted 
value from TW96 for the critical velocity for suppressing gas 
infall, 75 km/s, refers to v\ji defined at twice the pressureless 
collapse time, 2f e . This particular choice can be motivated by 
the time-scale of mergers between halos. Typical halos merge 
into more massive halos on a timescale that is approximately 



the Hubble time (Lacey & Cole 1993), so that requiring gas to 
cool and collapse by 2f e is roughly equivalent to requiring this 
to occur before the halo undergoes significant mergers (at which 
stage, gas can, in any case, cool and condense more effectively 
into the deeper potential well of the the newly formed merger 
product). In the context of reionization, we are interested in po- 
tentially more prompt feedback, since the global neutral frac- 
tion at high redshift, when the rare, high-sigma peaks are col- 
lapsing, can be evolving on a shorter time-scale ( <S 10% of the 
Hubble time, e.g., Haiman & Loeb 1997, 1998). It is therefore 
useful to know the amount of the baryons available immedi- 
ately after the collapse of the halo for star-formation, and thus 
for contributing to reionization. For completeness, the values 
of vi/2 and vo after t = t c are shown in Table 1, and compared to 
these values measured after 2f c . It can be seen from this table 
that the final results for the velocities (bottom row) remain un- 
changed. However, the relative importance of the three causes 
(§ 3) determining their value, and driving them to be smaller 
than in our z = 2 runs, are different (cooling contributing more 
when t c is used). 

4.4. The Impact on the Global Reionization History 

The presence or absence of a strong UV feedback on halos 
with circular velocities just above the threshold v c j rc <; lOkms -1 
has a strong influence on the global reionization history. As 
argued by Oh & Haiman (2003; see also Ricotti et al. 2002; 
Haiman, Abel & Rees 2000; Ciardi et al. 2000), halos with 
Vtirc lOkms -1 are unlikely to significantly contribute to reion- 
ization. As a result, the onset of the reionization epoch, when 
a significant fraction of the IGM can first be ionized, hinges 
on the efficiency of ionizing photon production in halos with 
Vdrc <J lOkms" 1 . Realistically, however, the UV feedback can- 
not have a global impact until a significant fraction of the IGM 
is already ionized (although see notes about clustering above). 
On the other hand, the UV feedback -if important - would im- 
pede the late stages of reionization and percolation, and delay 
the epoch when the Gunn-Peterson 'breakthrough' occurs. 

To illustrate this more clearly, we can make the following 
simple comparison. Let us assume that halos with circular ve- 
locities up to v c i rc = 75kms _1 are excluded from by strong feed- 
back from completing the reionization epoch (e.g., Haiman & 
Holder 2003; Wyithe & Loeb 2003). In this case, percolation 
will be achieved by halos with v c j rc > 75kms _1 . For example, 
using the mass function in Jenkins et al. (2001), we find that the 
universal fraction of baryons collapsed into halos of this size at 
z = 7 is ~ 1%. With reasonable choices of the efficiency pa- 
rameters (e.g. Haiman & Holder 2003; Cen 2003), the IGM 
will be fully ionized close to this redshift at 6 <J z ^ 7. In 
comparison, if we assume that all halos with v c j rc > lOkms" 1 
can contribute to reionization (representing the other extreme 
case of no UV feedback at all), then we find the collapsed frac- 
tion at z = 7 is ~ 10%, and the fraction of ~ 1% is achieved 
already by z = 12, making percolation feasible at this earlier 
epoch. In other words, under the assumption of fixed efficien- 
cies, the presence/absence of the feedback would make a fac- 
tor of ~ 1 .7 difference for the redshift of the "break-through" 
epoch. Assuming the IGM is suddenly fully ionized at this 
break-through, this translates to a factor of ~ 2.2 increase in 
the electron scattering optical depth, from r ~ 0.05 at z = 7 to 
t ~ 0.11 atz= 12. 

We found in this paper that UV feedback is unimportant at 
high redshift. Naively, this would suggest that relatively smaller 
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halos can contribute to reionization, and help explain the large 
optical depth measured by WMAP. On the other hand, as ex- 
plained in the introduction, the WMAP results suggest that the 
reionization history is more complex - with a "percolation" be- 
ing completed at redshifts of z ~ 6— 7, but with a significant tail 
of (perhaps only partial) ionization extending to much higher 
redshift. As clear from the preceding paragraph, a strong UV 
feedback would, in fact, have naturally helped to produce this 
behavior (as argued, e.g., in Haiman & Holder - see Figure 
4 in that paper). Indeed, provided halos near the <~ 10 km/s 
threshold otherwise have a high efficiency of ionizing photon 
production, a strong UV feedback would result in an extended 
episode during which the IGM is kept partially ionized. Full 
percolation could then be achieved only around z ~ 7 when ha- 
los with v c ; rc ~ 75kms _1 , impervious to feedback, start appear- 
ing. The absence of the UV feedback means that some other 
physics must then be responsible for the effective decrease in 
the global emissivity of ionizing radiation. Such a decrease is 
required to explain the complex reionization history suggested 
by the combination of WMAP and lower redshift data. 

5. CONCLUSIONS 

In this paper, we investigated the importance of photoioniza- 
tion heating in suppressing gas cooling and infall in halos with 
circular velocities in the range 10-100 km/s, collapsing at high 
redshift. We used a modified version of the one-dimensional, 
spherically symmetric, Lagrangian hydro code written by Thoul 
& Weinberg (1995), thereby directly extending the earlier cal- 
culations by Thoul & Weinberg (1996) to redshifts z > 10. We 
reproduced the results of TW96, namely that UV background 
has a strong impact in the low-redshift universe (z <J 3), sup- 
pressing gas infall by a factor of ~ 2 in galactic halos with cir- 
cular velocities of V1/2 ~ 50kms _1 . The main new result of 
the present paper is that suppression is much less significant at 
high redshift. For halos collapsing at z = 1 1 , we found vi /2 = 20 
km/s, and that a few percent of the gas can cool and condense 
in halos with circular velocities as low as vo=10 km/s. 

In § 4, we identified three reasons for this result, namely, 
(1) the amplitude of the ionizing background at high redshift 
is lower, (2) the ionizing radiation turns on when the perturba- 
tion that will become the dwarf galaxy has already grown to a 
substantial overdensity, and (3) collisional cooling processes at 
high redshift are more efficient. We find that the relative impor- 
tance of these three effects depends on the details of how and at 
what epoch the collapsed gas is identified, but we find that the 
combined effect is to essentially eliminate strong UV feedback 
in all halos with v c ; rc ~ lOkms -1 . 

In our simulations, the halos were assumed to be optically 
thin, and spherically symmetric. We explicitly showed that 
the halos are, in fact, self-shielding, and the radiation field is 
expected to be very anisotropic (the halos being illuminated 
by a single nearby source). Issues of radiative transfer and 
three-dimensional topology are therefore crucial in more de- 
tailed modeling of the UV feedback studied here. However, 
both effects will significantly strengthen our main result, namely 
that feedback is largely eliminated in high-redshift halos. 

The presence or absence of a strong photoheating feedback, 
which determines whether stars can or cannot form in halos just 
above the threshold v c j rc lOkms" 1 , has a strong influence on 
this reionization history. While the lack of a UV feedback, as 
suggested here, can help explaining the high electron scattering 
optical depth detected by WMAP, it leaves the puzzle that some 



other physics must be responsible for the complex reionization 
history suggested by combining the WMAP results with data 
from lower redshifts. 
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FIG. 1 . — The fraction / co u of gas that has collapsed as a function of circular velocity, normalized to the mass M c (p = 0) of gas that would have collapsed in the 
same time in the absence of any pressure. Upper left panel: Circles show halos collapsing at z = 2, while triangles represent halos that collapse at Zc = 1 1 ■ Feedback 
becomes much less important at higher redshifts. This is even more so when self shielding is considered: The triangle at v c j rc = 15 km/s moves from f co n = upto 
~ 1 when radiative transfer is taken into account. In the next three panels, we investigate the reasons for this difference, changing the parameters of the low-redshift 
model, one-by-one in succession, until we arrive back at the high redshift model. Upper right panel: increased cooling rate at higher redshift and densities. Lower 
left panel: lowering the amplitude of the UV-flux; lower right panel: late turn-on of the flux (see text for discussion). 
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FIG. 2. — The figure shows (upper to lower panels, respectively) the evolution of the radii, temperatures, ionization fractions, and cooling times, of shells during 
the collapse of a halo with Voire = 30 km/s. The set of panels on the left/right assumes Zon = 17/oo. Top row: The left panel shows radii of gas mass shells initially 
enclosing masses of 0.5Mf and Mf (shell 0.5 and 1), with J21 = 10~ 2 and in the absence of any pressure (pair of solid curves). The thick dotted curves illustrate 
the effect of a late turn-on for the flux. In the right panel, the dotted/dashed curves show the case with J 2 i = 10~ 2 /1. The solid curves, that show the pressureless 
collapse, are reproduced from the left panel. The panels below this all focus on shell 1, with the set on the left/right side showing the 721 = 10 _2 /1 case. Second row: 
The temperature evolution of the gas shell is shown together with the CMB temperature (solid curve), and an adiabatic temperature evolution (dot-dashed curve). 
Third row: The abundances of neutral and ionized hydrogen are plotted. Bottom row: Cooling timescales for the three dominant cooling processes; Compton cooling 
(solid curve), recombination (dashed curve) and excitation cooling (dotted curve). In the left panel, collisional cooling of neutral hydrogen dominates. Before z on 
the cooling time for this process is 2> 10 10 yrs and rapidly drops to a few times 10 7 beyond z on . In the right panel, Compton cooling dominates. Collisional cooling 
is very inefficient and is off the plot. 
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FIG. 3. — The evolution of V\n as a function of redshift. Circles denote the case with 721 = 10 2 for all redshifts greater than ion- Triangles denote the case when 
721 evolves as a function of redshift. 
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FIG. 4. — The optical depth from the virial radius inward to half this radius for 15 km/s NFW DM halos with their gas in hydrostatic equilibrium. The two curves 
denote the optical depth when the gas is in ionization equilibrium with an ionizing background 7„ = /2l(~?0 x 10~ 21 erg/s/cm 2 /sr/Hz. Two different values of 72i in 
the figure are shown, as labeled. The open dashed circles show the optical depth as extracted from the fiducial high redshift run. The optical depth is evaluated at 
the instant the ionizing background is switched on. 



